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A SYSTEM AND PROCESS FOR AUTOMATIC EXPOSURE CORRECTION IN 

AN IMAGE 

BACKGROUND 

Technical Field: 

The invention is related to improving the appearance of digitized or digital 
images, and more particularly to a system and process for improving the 
appearance of improperly exposed and/or improperly colored images. 

Background Art: 

Color images captured with a digital camera or obtained by scanning (i.e., 
digitizing) photographs can exhibit a less than appealing appearance owing to 
under or over exposed pixels and improperly colored pixels. In the case of 
exposure, capturing an image or taking a photograph using the proper exposure 
is often a difficult task. Several factors can contribute to poorly exposed images 
and photographs. For example, modern cameras often have automatic systems 
for estimating the proper exposure. However, these systems are sometimes 
fooled owing to the lighting conditions or atypical characteristics of the scene 
being captured. In addition, these modern cameras, as well as their 
predecessors, often allow the photographer to manually adjust the exposure 
setting, thereby introducing the possibility of human error. 

Improperly colored pixels typically occur because the colors of images 
recorded by a camera depend on the spectral characteristics of the camera 
optics, the lighting conditions of the scene and the characteristics of the color 



filters employed in the camera. As a result, the primary color data captured may 
not be a faithful rendition of the actual colors found in the scene. 

To resolve the problems of improperly exposed and colored images, post- 
5 capture corrective measures are often employed. There are two categories of 
systems. In the first category, the systems that implement these corrective 
measures can be manual in that, for example, a user views the image on a 
computer screen and employs a user interface to interactively adjust the color 
and/or exposure. In the second category, the systems provide for automatic 

10 correction of color and/or exposure in an image. In these systems, a user 

typically just accesses the image, which can be viewed on a computer screen, 
and activates an automatic correction program. This activation can be implicit, in 
the sense that the automatic correction program is called whenever the image is 
accessed; it can also be explicit, in which case the user clicks on a button to call 

15 the correction program when necessary. The present invention falls into this 
latter automatic correction category. 

SUMMARY 

20 

The present invention is directed toward a system and process for 
improving the appearance of improperly exposed and/or improperly colored 
images. This involves the use of two novel techniques-namely an automatic 
color correction technique and an automatic exposure correction technique. 
25 The automatic color correction technique takes information from within an image 
to determine true color characteristics, and improves the color in improperly 
colored pixels. The automatic exposure correction technique measures the 
average intensity of all of the pixels and adjusts the entire image pixel by pixel to 
compensate for over or under exposure. 
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In regard to the automatic color correction technique, this generally 
involves a system and process for correcting the color of improperly colored 
pixels of an image having multi-level primary color channels (e.g., red (R), green 
(G), and blue (B)) by first computing a histogram of the color levels of the image 
5 pixels for each of the color channels. The lowermost and uppermost color levels 
for each of the color channels, which are consistent with the overall distribution 
of color levels for the corresponding channel, are then computed using the color 
channel histograms. In addition, the per channel average color level for a group 
of the brightest pixels is computed. The lowermost and uppermost color levels 

10 for each of the color channels and the per channel average color levels of the 
brightest pixels are used to compare the color levels of pairs of the color 
channels to determine if the color levels in each of the pairs is balanced. In 
tested embodiments of the present automatic color correction system and 
process employing RGB color channels, the green color channel is used as the 

15 reference for balancing such that the green and blue color channels are 

compared and the green and red color channels are compared. This was done 
because in the Bayer CCD chips typically used in digital cameras there are twice 
as many G pixels than R and B pixels, and so the G values are considered to be 
more reliable. However, this assumption need not be made and another of the 

20 color channels can be used as the reference channel. A pair of color channels is 
considered balanced if the dynamic range of their color levels match within a 
prescribed tolerance. If the color channels are balanced, then no correction is 
required. However, if a compared channel pair is found to be out of balance, 
then the dynamic range of the color channel with the narrower range is linearly 

2 5 expanded to match the channel with the wider dynamic range to a desired 
degree. 

In regard to ascertaining if the dynamic color level ranges of a compared 
channel pair match within a prescribed tolerance, this involves first determining if 
30 the difference between the average color levels of the group of the brightest 
pixels for the channel pair under consideration is less than their average 
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multiplied by a prescribed tolerance value. It is also determined whether the 
difference between the lowermost color levels of color channel pair is less than 
their average multiplied by the prescribed tolerance value. If either of the 
computed differences is determined not to be less, the dynamic color level 
5 ranges of the channel pair being compared is deemed not to match. 

The aforementioned desired degree of matching between the dynamic 
range of the color channel with the narrower range and the range of the channel 
with the wider range refers to a potential tailoring of how aggressively the color 

10 correction is performed. More particularly, this involves varying the degree to 
which the color channel with the narrower range is made to match the channel 
with the wider dynamic range between the extremes of expanding the channel 
with the narrower range to match the full range of the channel with the wider 
range, to no expansion at all, depending on the value of a safeness parameter. 

15 This safeness parameter can be prescribed or user-specified. 

If the safeness parameter is employed, linearly expanding the dynamic 
range of the color channel with the narrower range to match the channel with the 
wider dynamic range to the desired degree involves first computing new 

2 0 lowermost and uppermost levels for the channel with the wider dynamic color 
level range that reflect the degree to which it is desired for the narrower range 
channel to match the wider range channel. The new lowermost and uppermost 
levels of the channel with the wider dynamic color level range and the lowermost 
and uppermost levels of the channel with the narrower range are then employed 

25 to establish a linear color correction transform that maps each old color level to a 
new color level. This linear color correction transform is then applied to the 
narrower range color channel to produce color corrected pixel levels for each 
pixel of the image. Due to use of the safeness parameter, the color correction 
program can be called multiple times to achieve more aggressive color 

30 correction. The color correction process will always converge. However, in one 
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embodiment, the safeness parameter was chosen such that one execution is 
enough to result in satisfactory color correction for a majority the tested images. 

In regard to the automatic exposure correction technique, this generally 
5 involves a system and process for correcting the exposure of improperly exposed 
pixels of an image by first linearly expanding the dynamic range of the intensity 
levels of the image pixels so as to match the full dynamic intensity range 
available for the application to a desired degree. It is then determined whether 
the linearly expanded intensity levels of the image pixels are evenly distributed. 
10 If they are, then no further correction is required. However, if the levels are not 
found to be evenly distributed, then a gamma correction factor is computed and 
applied to the linearly expanded intensity level of each pixel to produce corrected 
intensity values. 

15 More particularly, the aforementioned linear expansion involves first 

computing a histogram of the intensity levels of the image pixels, and then 
computing the lowermost and uppermost intensity levels that are consistent with 
the overall distribution of intensity levels in the image. Next, new lowermost and 
uppermost intensity levels are computed that span the full dynamic intensity 

20 range available to the desired degree. The originally computed lowermost and 
uppermost intensity levels and the new lowermost and uppermost intensity levels 
are employed to establish a linear intensity correction transform that maps the 
original intensity level of each pixel in the image to a linearly expanded intensity 
level. The linear intensity correction transform is then to each pixel of the image. 

25 

The aforementioned desired degree of matching between the dynamic 
range of the original intensity levels and the full dynamic intensity range available 
refers to a potential tailoring of how aggressively the exposure correction is 
performed. More particularly, this involves varying the degree to which intensity 
30 levels of the image pixels are made to match the full dynamic intensity range 

available between the extremes of expanding the dynamic range of the intensity 
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levels to match the full dynamic intensity range, to no expansion at all, 
depending on the value of a safeness parameter. This safeness parameter can 
be prescribed or user-specified, and need not be the same parameter used with 
the previously-described color correction technique. Due to use of the safeness 
parameter, the exposure correction program can be called multiple times to 
achieve more aggressive exposure correction. The exposure correction process 
will always converge. However, in one embodiment, the safeness parameter 
was chosen such that one execution is enough to produce satisfactory exposure 
correction for the majority of the tested images. 

In regard to determining whether the linearly expanded intensity levels of 
the image pixels are evenly distributed, this generally involves first computing a 
histogram from the linearly expanded pixel intensity levels. A 50 percentile (i.e., 
median) intensity level is then computed from the linearly expanded pixel 
intensity level histogram as the level where (1) the sum of the number of pixels 
exhibiting a level within a prescribed range of levels extending from lowest level 
possible to the unknown median level, less one level, is less than the total 
number of pixels in the image multiplied by 0.5, and (2) the sum of the number of 
pixels exhibiting a level within a prescribed range of levels extending from lowest 
level possible to the unknown median level, is greater than or equal to the total 
number of pixels in the image multiplied by 0.5. It is next determined if the 
difference between the median intensity level and one half the maximum 
intensity level is less than or equal to the maximum intensity value multiplied by 
an intensity tolerance factor. If it is, then the intensity levels are evenly balanced 
and no further exposure correction is required. Otherwise, the linearly expanded 
intensity levels are deemed not to be evenly distributed. 

In cases where the linearly expanded intensity levels are deemed not to 
be evenly distributed, a nonlinear gamma correction factor is computed as 
indicted previously. This generally involves first determining whether there is a 
significant peak in the lower side, upper side, or both sides of the linearly 



expanded pixel intensity level histogram. A peak is considered significant if the 
number of pixels assigned to the bins of the histogram around the peak pixel 
intensity level is larger than the average. A significant peak determined to be in 
the lower side of the histogram, and which exhibits the lowest intensity level in 
5 cases where there are more than one, is designated as the low peak. Similarly, 
a significant peak determined to be in the upper side of the histogram, and which 
exhibits the highest intensity level in cases where there are more than one, is 
designated as the high peak. In addition to identifying the low and high peaks, 
the 25 percentile and 75 percentile intensity levels are computed from the 

10 linearly expanded pixel intensity level histogram, and used along with the 

previously computed 50 percentile (i.e. median) intensity level, to establish the 
gamma correction factor. More particularly, the gamma correction factor is 
established using the foregoing information as follows. When the median 
intensity level is greater than the maximum possible level multiplied by 0.50 plus 

1 5 the intensity tolerance factor: 

(1) the correction factor is set to 1 if the low peak is less than one- 
eighth the maximum possible intensity level and the median intensity level is 
greater than or equal to three-quarters the maximum possible intensity level; 

(2) the correction factor is set to the log of 0.125 divided by the log of 
20 the low peak divided by the maximum possible intensity level if the low peak is 

less than one-eighth the maximum possible intensity level and the median 
intensity level is less than three-quarters the maximum possible intensity level; 

(3) when the low peak is greater than one-eighth the maximum 
possible intensity level and the median intensity level is greater than three- 

25 quarters the maximum possible intensity level, the correction factor is 

interpolated between 1 and a prescribed maximum value if the 75 percentile 
intensity level is between the maximum possible intensity level, and seven- 
eighths the maximum possible intensity level; 

(4) when the low peak is greater than one-eighth the maximum 
30 possible intensity level and the median intensity level is greater than three- 
quarters the maximum possible intensity level, and 75 percentile intensity level is 
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less than or equal to seven-eighths the maximum possible intensity level, the 
correction factor is either set to the log of 0.5 divided by the log of the median 
intensity level divided by the maximum possible intensity level, or in an alternate 
embodiment set to the log of 0.55 divided by the log of the median intensity level 
divided by the maximum possible intensity level except that the factor is limited 
to at most the prescribed maximum correction value. This alternate embodiment 
generally reduces the difference between the corrected image and the original 
image, thereby providing a more psychologically pleasing correction to a typical 
user; and 

(5) in all other cases, the correction factor is set to 1 to indicate that no 
correction is necessary. 

Conversely, when the median intensity level is less than the maximum 
possible level multiplied by 0.50 less the intensity tolerance factor: 

(1 ) the correction factor is set to 1 if the high peak is greater than 
seven-eighths the maximum possible intensity level; 

(2) the correction factor is set to 7 = 1 - 0 - 7 min ) whenever the 25 

percentile intensity level is less than or equal to one-eighth the maximum 
possible intensity level, where y is the correction factor, L is the maximum 
possible intensity level, Q0.25 is the 25 percentile intensity level and y min is a 
prescribed minimum correction factor value; and 

(3) in all other cases, the correction factor is set to the log of 0.5 
divided by the log of the median intensity level divided by the maximum possible 
intensity level, except that the factor is limited to at least the prescribed minimum 
correction factor value. 

It is noted that whenever gamma correction results in an overall 
brightening of the image, the color saturation levels of each color channel of 
each pixel of the image can be boosted by a desired amount to offset a typically 
perceived effect that the colors have been washed out. This boosting can entail 
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varying the degree to which the color saturation levels of the image pixels are 
boosted between about 25 percent to less than 1 percent, in proportion to the 
degree of brightening. 

It is further noted that while tested embodiments of present exposure 
correction system and process employed the foregoing gamma function, other 
functions such as splines, can be used as well without significant adverse effect. 

In addition to the just described benefits, other advantages of the present 
invention will become apparent from the detailed description which follows 
hereinafter when taken in conjunction with the drawing figures which accompany 
it. 

DESCRIPTION OF THE DRAWINGS 

The specific features, aspects, and advantages of the present invention 
will become better understood with regard to the following description, appended 
claims, and accompanying drawings where: 

FIG. 1 is a diagram depicting a general purpose computing device 
constituting an exemplary system for implementing the present invention. 

FIGS. 2(a)-(b) are images illustrating the results of the application of the 
present automatic exposure correction technique. Fig. 2(a) depicts a seagull in 
an image that is underexposed. Fig. 2(b) shows the same image after having 
the exposure corrected. 

FIG. 3 is a flow chart diagramming an overall process for correcting the 
exposure of an image in accordance with the present invention. 



FIGS. 4(a)-(b) are images illustrating the results of the application of the 
present automatic color correction technique. Fig. 4(a) depicts a baby in an 
image that is improperly colored. Fig. 4(b) shows the same image after the color 
is corrected. 

FIG. 5 is a flow chart diagramming an overall process for correcting the 
color of an image in accordance with the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the following description of the preferred embodiments of the present 
invention, reference is made to the accompanying drawings which form a part 
hereof, and in which is shown by way of illustration specific embodiments in 
which the invention may be practiced. It is understood that other embodiments 
may be utilized and structural changes may be made without departing from the 
scope of the present invention. 

1.0 The Computing Environment 

Before providing a description of the preferred embodiments of the 
present invention, a brief, general description of a suitable computing 
environment in which the invention may be implemented will be described. 
Figure 1 illustrates an example of a suitable computing system environment 100. 
The computing system environment 100 is only one example of a suitable 
computing environment and is not intended to suggest any limitation as to the 
scope of use or functionality of the invention. Neither should the computing 
environment 100 be interpreted as having any dependency or requirement 
relating to any one or combination of components illustrated in the exemplary 
operating environment 100. 
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The invention is operational with numerous other general purpose or 
special purpose computing system environments or configurations. Examples of 
well known computing systems, environments, and/or configurations that may be 
suitable for use with the invention include, but are not limited to, personal 
computers, server computers, hand-held or laptop devices, multiprocessor 
systems, microprocessor-based systems, set top boxes, programmable 
consumer electronics, network PCs, minicomputers, mainframe computers, 
distributed computing environments that include any of the above systems or 
devices, and the like. 

The invention may be described in the general context of computer- 
executable instructions, such as program modules, being executed by a 
computer. Generally, program modules include routines, programs, objects, 
components, data structures, etc. that perform particular tasks or implement 
particular abstract data types. The invention may also be practiced in distributed 
computing environments where tasks are performed by remote processing 
devices that are linked through a communications network. In a distributed 
computing environment, program modules may be located in both local and 
remote computer storage media including memory storage devices. 

With reference to Figure 1 , an exemplary system for implementing the 
invention includes a general purpose computing device in the form of a computer 
110. Components of computer 110 may include, but are not limited to, a 
processing unit 120, a system memory 130, and a system bus 121 that couples 
various system components including the system memory to the processing unit 
120. The system bus 121 may be any of several types of bus structures 
including a memory bus or memory controller, a peripheral bus, and a local bus 
using any of a variety of bus architectures. By way of example, and not 
limitation, such architectures include Industry Standard Architecture (ISA) bus, 
Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video 



Electronics Standards Association (VESA) local bus, and Peripheral Component 
Interconnect (PCI) bus also known as Mezzanine bus. 

Computer 1 1 0 typically includes a variety of computer readable media. 
Computer readable media can be any available media that can be accessed by 
computer 1 10 and includes both volatile and nonvolatile media, removable and 
non-removable media. By way of example, and not limitation, computer 
readable media may comprise computer storage media and communication 
media. Computer storage media includes both volatile and nonvolatile, 
removable and non-removable media implemented in any method or technology 
for storage of information such as computer readable instructions, data 
structures, program modules or other data. Computer storage media includes, 
but is not limited to, RAM, ROM, EEPROM, flash memory or other memory 
technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, 
magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic 
storage devices, or any other medium which can be used to store the desired 
information and which can be accessed by computer 110. Communication 
media typically embodies computer readable instructions, data structures, 
program modules or other data in a modulated data signal such as a carrier 
wave or other transport mechanism and includes any information delivery media. 
The term "modulated data signal" means a signal that has one or more of its 
characteristics set or changed in such a manner as to encode information in the 
signal. By way of example, and not limitation, communication media includes 
wired media such as a wired network or direct-wired connection, and wireless 
media such as acoustic, RF, infrared and other wireless media. Combinations of 
the any of the above should also be included within the scope of computer 
readable media. 

The system memory 130 includes computer storage media in the form of 
volatile and/or nonvolatile memory such as read only memory (ROM) 131 and 
random access memory (RAM) 132. A basic input/output system 133 (BIOS), 
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containing the basic routines that help to transfer information between elements 
within computer 110, such as during start-up, is typically stored in ROM 131. 
RAM 132 typically contains data and/or program modules that are immediately 
accessible to and/or presently being operated on by processing unit 120. By way 
of example, and not limitation, Figure 1 illustrates operating system 134, 
application programs 135, other program modules 136, and program data 137. 

The computer 110 may also include other removable/non-removable, 
volatile/nonvolatile computer storage media. By way of example only, Figure 1 
illustrates a hard disk drive 141 that reads from or writes to non-removable, 
nonvolatile magnetic media, a magnetic disk drive 151 that reads from or writes 
to a removable, nonvolatile magnetic disk 152, and an optical disk drive 155 that 
reads from or writes to a removable, nonvolatile optical disk 156 such as a CD 
ROM or other optical media. Other removable/non-removable, 
volatile/nonvolatile computer storage media that can be used in the exemplary 
operating environment include, but are not limited to, magnetic tape cassettes, 
flash memory cards, digital versatile disks, digital video tape, solid state RAM, 
solid state ROM, and the like. The hard disk drive 141 is typically connected to 
the system bus 121 through an non-removable memory interface such as 
interface 140, and magnetic disk drive 151 and optical disk drive 155 are 
typically connected to the system bus 121 by a removable memory interface, 
such as interface 150. 

The drives and their associated computer storage media discussed above 
and illustrated in Figure 1 , provide storage of computer readable instructions, 
data structures, program modules and other data for the computer 110. In 
Figure 1, for example, hard disk drive 141 is illustrated as storing operating 
system 144, application programs 145, other program modules 146, and 
program data 147. Note that these components can either be the same as or 
different from operating system 134, application programs 135, other program 
modules 136, and program data 137. Operating system 144, application 
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programs 145, other program modules 146, and program data 147 are given 
different numbers here to illustrate that, at a minimum, they are different copies. 
A user may enter commands and information into the computer 110 through 
input devices such as a keyboard 162 and pointing device 161, commonly 
5 referred to as a mouse, trackball or touch pad. Other input devices (not shown) 
may include a microphone, joystick, game pad, satellite dish, scanner, or the like. 
These and other input devices are often connected to the processing unit 120 
through a user input interface 160 that is coupled to the system bus 121, but 
may be connected by other interface and bus structures, such as a parallel port, 

10 game port or a universal serial bus (USB). A monitor 191 or other type of display 
device is also connected to the system bus 121 via an interface, such as a video 
interface 190. In addition to the monitor, computers may also include other 
peripheral output devices such as speakers 197 and printer 196, which may be 
connected through an output peripheral interface 195. Of particular significance 

15 to the present invention, a camera 1 92 (such as a digital/electronic still or video 
camera, or film/photographic scanner) capable of capturing a sequence of 
images 193 can also be included as an input device to the personal computer 
110. Further, while just one camera is depicted, multiple cameras could be 
included as input devices to the personal computer 110. The images 193 from 

20 the one or more cameras are input into the computer 1 10 via an appropriate 
camera interface 194. This interface 194 is connected to the system bus 121, 
thereby allowing the images to be routed to and stored in the RAM 132, or one of 
the other data storage devices associated with the computer 110. However, it is 
noted that image data can be input into the computer 1 10 from any of the 

25 aforementioned computer-readable media as well, without requiring the use of 
the camera 192. 

The computer 110 may operate in a networked environment using logical 
connections to one or more remote computers, such as a remote computer 180. 
30 The remote computer 180 may be a personal computer, a server, a router, a 

network PC, a peer device or other common network node, and typically includes 
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many or all of the elements described above relative to the computer 110, 
although only a memory storage device 181 has been illustrated in Figure 1. 
The logical connections depicted in Figure 1 include a local area network (LAN) 
171 and a wide area network (WAN) 173, but may also include other networks. 
5 Such networking environments are commonplace in offices, enterprise-wide 
computer networks, intranets and the Internet. 

When used in a LAN networking environment, the computer 110 is 
connected to the LAN 171 through a network interface or adapter 170. When 

10 used in a WAN networking environment, the computer 1 10 typically includes a 
modem 172 or other means for establishing communications over the WAN 173, 
such as the Internet. The modem 172, which may be internal or external, may 
be connected to the system bus 121 via the user input interface 160, or other 
appropriate mechanism. In a networked environment, program modules 

15 depicted relative to the computer 1 10, or portions thereof, may be stored in the 
remote memory storage device. By way of example, and not limitation, Figure 1 
illustrates remote application programs 185 as residing on memory device 181. 
It will be appreciated that the network connections shown are exemplary and 
other means of establishing a communications link between the computers may 

20 be used. 

2.0 The Automatic Color And Exposure Correction System And Process 

The exemplary operating environment having now been discussed, the 
25 remaining part of this description section will be devoted to a description of the 
program modules embodying the invention. Generally, the system and process 
according to the present invention involves two techniques for users to improve 
the appearance of imperfectly taken digital images or digitized photographs or 
frames of a video (hereinafter collectively referred to as an image). These 
30 techniques are automatic exposure correction and automatic color correction. 
The automatic exposure correction technique generally measures the average 
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intensity of all of the pixels and can adjust the entire image pixel by pixel to 
compensate for over- or under-exposure. The automatic color correction 
technique generally takes information from within an image to determine true 
color characteristics, and improves the color in improperly colored pixels. These 
5 techniques are stand alone in that each can be applied to an image exclusive of 
the other, although it is envisioned that they would be packaged together such 
that a user activates one program which then automatically corrects both the 
exposure and color of an image. It does not matter which technique is applied 
first, however it is noted that in tested embodiments color correction was done 
10 first followed by exposure correction. 

2.1 Automatic Exposure Correction 

Automatic exposure correction according to the present system and 
15 process is used to better reveal the details captured in an image which is under 
or over exposed. The technique is based on histogram processing, and 
generally involves an initial linear transformation, which is followed, if necessary, 
with a nonlinear gamma correction in order to achieve the best result. In some 
cases a color saturation boost procedure is also used. The results of the 
20 application of the present automatic exposure correction technique is shown in 
the before and after images found in Figs. 2(a)-(b). Fig. 2(a) depicts a seagull in 
an image that is underexposed. Fig. 2(b) shows the same image after having 
the exposure corrected in accordance with the present technique. 

25 More particularly, as outlined in the flow diagram in Fig. 3, the general 

process for implementing automatic exposure correction in accordance with the 
present invention involves obtaining a digitized or digital image which a user 
wants to improve by correcting the exposure levels of the image pixels (process 
action 300). The dynamic range of the intensity levels of the image pixels is then 

30 linearly expanded so as to match the full dynamic intensity range available to a 
desired degree (process action 302). The desired degree of matching refers to 
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potential tailoring of how aggressively the exposure correction is performed and 
will be discussed in more detail shortly. Once the pixel intensity levels have 
been linearly expanded, it is determined if the expanded levels are evenly 
distributed (process action 304). If so, no further correction is required and the 
5 process ends. However, if it is found that the levels are not evenly distributed, 
then a non-linear gamma correction is performed by computing a gamma 
correction factor (process action 306) and applying it to each linearly expanded 
pixel intensity level (process action 308) to further correct the exposure of each 
pixel of the image. Optionally, once the gamma correction is complete, it can be 

10 determined if the exposure correction has resulted in an overall brightening of 
the image (optional process action 310). If not, the process ends. However, if 
the image was brightened, the color saturation levels of each color channel of 
each pixel of the image is boosted (optional process action 312). In tested 
embodiments of the exposure correction procedure, the degree of boosting was 

15 varied depending on the amount of brightening as will also be explained in more 
detail below. It is noted that the optional nature of the last two process actions is 
indicated by using dashed lines in the flow diagram of Fig. 3. 

Each of the foregoing phases of the present exposure correction system 
20 and process will now be described separately in the sub-sections to follow. 

2.1.1 Linear Intensity Correction 

A histogram of the pixel intensities is computed as a discrete function h(r k ) 
25 = n k , where r k is the Mh intensity level and n k is the number of pixels in the 

image having intensity level r k . Given the intensity histogram, the lower (bottom) 
and upper (top) intensity levels, denoted by B and T, respectively, are computed. 
It is noted these are not the actual highest and lowest values as outlying intensity 
values may exist in the image being corrected. Rather, the lower and upper 
30 intensity levels represent the lowest and highest intensity levels that are still 
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consistent with the overall distribution of intensity values in the image. More 
particularly, B is determined such that, 

fl-l B 

Y,n,<pN and Y i n l >pN i M) 

/=0 /=0 

where / is the intensity level in a range [0, L), N is the total number of pixels, and 
p is a small number chosen to eliminate the impact of noise at the top and 
bottom of the intensity range including eliminating saturated pixels. Similarly, T 
is determined such that, 

2_,n,>pN and < pN ' . (2) 

In tested embodiments of the present system and process employing a pixel 
intensity range of 0 to 255 integer levels, p was empirically chosen to be about 
0.5%. This represents an amount just slightly larger than one intensity level. It is 
noted that if more than 256 intensity levels are employed, the value of p should 
remain about 0.5%. However, if less than 256 intensity levels are employed, 
then the value of p should be adjusted to represent about 1 intensity level. For 
example, if 64 intensity levels are employed, the value of p should be about 
0.2%. 

In designing the present system and process, it was assumed that a well- 
exposed image is defined as one that, among other things, uses the whole 
dynamic range [0, L]. Thus, once the existing values of the lowest and highest 
intensity level (i.e., B and T) are computed, the expected values for these levels 
after correction (denoted by B and f), is computed assuming the intensity 
range is expanded. This is accomplished using the following equations: 
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B = (\-a)B and 



(3) 



f = a(L-l) + (l-a)T (4) 

The factor a represents a so-called safeness parameter that controls how 
conservatively the correction is performed. No correction will be performed if a = 
0, and the maximum correction will be performed if a = 1 . This factor can be 
prescribed, or user specified. If prescribed, since this is an automatic tool, the 
value should be chosen so as to be conservative in the correction. This may not 
achieve the best correction, but should have the effect of not making any image 
look worse. For example, a prescribed value for a of about 0.75 would be a 
good conservative choice. However, a user specified value may be a better way 
go since there will be exceptions to the assumption that an image should fill the 
entire dynamic intensity range available because the photographer may want to 
achieve certain effects. For example, in a portrait it is common to make a large 
fraction of the background black. Thus, expanding the dynamic range may 
lighten the background unacceptably. In such a case, the user would specify a 
low value fora . Thus, a user specified safeness parameter provides a way to 
tailor the amount of correction to the circumstances of the image. 

Once the existing and expected values for the high and low intensity 
levels have been computed, a linear intensity correction transformation is 
established and each old intensity level / is mapped to a new intensity level / . 
Specifically, the transformation takes the form of: 



l=aUb t (5) 



where, 



a = (f-B)/(T-B), and 



(6) 



b = (BT-TB)/(T-B) (?) 

2.1.2 Gamma Correction 

Another assumption made in designing the present system and process is 
that a well-exposed image is one in which the bright and dark pixels are 
distributed evenly; that is, the median of the intensity in a picture is expected to 
be around L/2. If this is not the case, then a correction is made. To accomplish 
this task it is first determined if the pixels are evenly distributed after completing 
the previously described linear transformation. If they are not, then a nonlinear 
gamma correction is applied in order to achieve the best overall exposure. 
Essentially, for generally overexposed images, the gamma correction has the 
effect of darkening the higher intensity pixel more than lower intensity pixels. 
This results in the brighter overexposed areas being significantly darken, while 
the darker areas are less affected so that there is no substantial loss of detail in 
the darker areas of the image. Conversely, for generally under exposed images, 
the gamma correction has the effect of brightening the lower intensity pixels 
more than the higher intensity pixels. This results in the darker underexposed 
areas being significantly brightened, while the brighter areas are less affected to 
avoid saturating the pixels in the bright areas. 

The distribution of pixels exhibiting darker and lighter intensities is 

A 

assessed as follows. First, a linearly transformed histogram K r k) is computed 
such that, 

h(ar k +b) = h(r k ) , Q) 
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Next, the 50-percentile (median) intensity level is computed, denoted by So.sfrom 
*( r *). More particularly Q p (p = 0.5) is determined such that, 

Q fkD< P N i3nd 

/=0 1 

5 f^KQZpN 

1=0 

It is then determined if, 

|Q 05 -£/2| <0.05Z 

If so, the distribution of darker and lighter pixels is deemed to be well balanced 
10 and processing stops. It is noted a factor of 0.05 was chosen empirically to 
provide a small tolerance range, which even given noise there is a high 
confidence that the darker and lighter pixels are balanced. As such, other values 
for this factor should be chosen if the noise characteristics of the image dictates 
a larger value or allow for a smaller value. It is also noted that the value of this 
15 factor is dependent of the number of intensity levels, which in the case of the 
tested embodiments was 256. If more or less intensity levels are employed, the 
value of the factor should be adjusted proportionally. 



(9) 
(10) 

(11) 



If, however, Eq. 1 1 is found not to be the case, the aforementioned 
20 nonlinear gamma correction is applied as follows. First, it is determined if there 
is a significant peak in the lower side of the linearly transformed histogram, or a 
significant peak in the upper side of this histogram. A peak is considered 
significant if the number of pixels around the peak is larger than the average. A 
peak P is significant if and only if, 

25 

2>(/)>— r-tf, (12 ) 
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where D is the neighborhood size. In operation, Eq. 12 is used to check each 
intensity level for a peak. The peak (if any) that is associated with the lowest 
intensity level is designated as P\ 0VI and the peak (if any) associated with the 
highest intensity level is designated as P h igh- 

5 

The aforementioned neighborhood size is chosen so as to be large 
enough exclude noise in any one intensity level from skewing the results, while at 
the same time not so large that a actual peak is lost in the averaging across 
intensity levels. In tested embodiments using 266 intensity levels, a value of D 
10 between 3 and 5 was found to be acceptable with 4 being preferred. It should be 
noted, however, that if fewer or more intensity levels are employed, the value of 
D should be adjusted proportionally. 

In addition to finding the peaks, the 25-percentile and 75-percentile 
15 intensity levels are computed, denoted by Qoisand Qo.is , respectively, from Kh) . 
More particularly Q p (p = 0.25 or 0.75) is determined using Eqs. (9) and (10). 
Thus, there is now a value computed for Qo.is , Qo.s and Qo.is . 

In cases where Qo.s > 0.551, an appropriate gamma y is established as 
20 follows, depending on the circumstances. 

1 . If the lower peak P iow < L/8 and Qo.s > 31/4, set y = 1 because it is hard to 
correct both ends with one gamma correction function. 

25 2. If the lower peak P iow < L/8, but Qo.s < 3L4, set y = log(l / 8) / log(/> ow / L) so 
as to map the lower peak to L/8 to make the details around the lower peak more 
visible, at the expense of making the overall image even brighter. 
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3. If the lower peak P l0VI 7z L/Q and Qo.$ > ZLIA, and also 60.75 > 7Z/8, the 
photographer may have some special intention to make the picture very bright. 

To account for this, y is interpolated between 1 (i.e., no correction) and Ymax as 
follows: 



r=l+Lr uT (r ™- l) - (13) 



Here, is the maximum value allowed in for this gamma correction, and can 

range between 1 .2 and 2.5. In tested embodiments Ymax was set to 1 .5 with 
10 satisfactory results. 

4. If the lower peak P\ ov/ > LI8 and Qo.$ > 3L/4, and also Qo.is < 7L/8, gamma 
can be set so that: 

15 r = log(O.5)/log(0 as /I) (14) 

so as to map Q0.5 to LI2. However, in one embodiment of the present system 
and process, gamma is set so that: 

20 y = log(0.55)/log(Q, 5 / L) (15) 

(note the "0.55" in the equation) to make the corrected image brighter than if the 
Eq. (14) is used. Furthermore, y is capped by/max . It was found that these 
modifications to Eq. (14) reduced the difference between the corrected image 
25 and the original image. Thus, while the gamma correction using Eq. (15) may 
not be optimal, the change it causes is less dramatic, thereby providing a more 
psychologically pleasing correction to the typical user. 
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4. If none of the foregoing special cases apply, set y = 1 . 



Conversely, in cases where Qo.s < 0.45L, an appropriate gamma y is 
established as follows, depending on the circumstances. 

1 . If the upper peak P histl > 71/8, set y = 1 because the upper peak would be 
washed out if an attempt is made to make the image brighter. 

2. If So w £ 1/8, the photographer may have some special intention to make 
the picture very dark. To account for this, y is set to 1 (no correction) if Qoas = 0 
and y is set to y mi „ if Qois = 1/8. For values of Qo is between 0 and 1/8, y is 
interpolated between Y mi n and 1. More particularly, the following equation can be 
used to compute y. 

Here, is the minimum value allowed in for this gamma correction, and can 

range between 0.4 and 0.56. In tested embodiments Ymin was set to 0.48 with 
satisfactory results. 

3. If none of the foregoing special cases apply, as in the case where Qo.s > 
0.551, gamma can be set so that: 

r = log(0.5)/log(& 5 /£) (17) 
to map 0o.5 to U2. However, if this results in r < , then r is set equal to Ymin. 
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Once the gamma correction factor r has been established, it is applied to 

correct each linearly transformed pixel intensity values / using the following 
equation: 



It is noted that while tested embodiments of present exposure correction 
system and process employed the foregoing gamma function, other functions 
such as splines, can be used as well without significant adverse effect. 



2.1.3 Color Saturation Boost 

Although the foregoing gamma correction does not alter the values of hue 
and saturation of the image, it usually affects the relative appearance of the 

15 colors (i.e., the overall color perception), especially when the image is made 
brighter. Specifically, the colors appear to wash out when the image is 
brightened. This perception issue can be resolved by boosting the color 
saturation levels of the pixels in an image whenever y < 1 , thereby indicating 
that the gamma correction brightened the image. In tested embodiments, it was 

2 0 found that an increase in color saturation levels of between 0 and 25 %, 

depending on the value of the gamma correction factor computed for the image 
produced the desired effect of deepening the appearance of the colors. 
Specifically, a full 25% boost is applied to each pixel if r = r mia t anc j a 
progressively smaller boost is applied the closer rgets to 1. In operation, this is 

2 5 accomplished by applying the following equation to each color channel of each 
pixel in the image: 



5 




(18) 



10 



1-y 



xl.25. 



(19) 



2.2 Automatic Color Correction 

Automatic color correction according to the present system and process 
5 takes information from within a photo to determine true color characteristics, and 
improves the color in improperly colored pixels. The technique is based on 
histogram processing for each of the primary color channels (e.g., R, G, B) and 
the luminance (Y) channel, and generally involves a linear transformation 
designed to balance the dynamic color range of the color channels. The results 
10 of the application of the present automatic color correction technique is shown in 
the before and after images found in Figs. 4(a)-(b). Fig. 4(a) depicts a baby in an 
image that is improperly colored. Fig. 4(b) shows the same image after the color 
is corrected in accordance with the present technique. 

15 More particularly, as outlined in the flow diagram in Fig. 5, the general 

process for implementing automatic color correction in accordance with the 
present invention involves obtaining a digitized or digital image which a user 
wants to improve by correcting the color levels of the image pixels (process 
action 500). A histogram of the color levels of the image pixels for each of the 

20 color channels is then computed (process action 502), followed by computing the 
lowermost and uppermost color levels for each of the color channels that are 
consistent with the overall distribution of color levels for that channel (process 
action 504). Next, the per channel average color level for a group of the 
brightest pixels is computed (process action 506). The lowermost and 

25 uppermost color levels for each of the color channels and the per channel 

average color levels are used to compare the color levels of pairs of the color 
channels (process action 508) to determine if they are balanced (process action 
510). A pair of color channels is considered balanced if the dynamic range of 
their color levels match within a prescribed tolerance, as will be described in 

30 more detail shortly. If the color channels are balanced, then no correction is 

required and the process ends. However, if any channel pair is found to be out 
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of balance, then the dynamic range of the color channel with the narrower range 
is linearly expanded to match the channel with the wider dynamic range to a 
desired degree (process action 512). The desired degree of matching refers to 
potential tailoring of how aggressively the color correction is performed and will 
be discussed in more detail below. 

Each of the foregoing phases of the present color correction system and 
process will now be described in more detail. 

A histogram for each of the primary color channels of the pixels in the 
image (which for this discussion will be assumed to be the R, G, B channels) is 
computed as discrete functions h R (r k ) = n k , h G {r k ) = n k , h B (r k ) = n ki where h R is the 
histogram for the red color channel, h G is the histogram for the green color 
channel, h B is the histogram for the blue color channel, r k is the k\h color level in 
one of the color channels and n k is the number of pixels in the image having 
color level r k for the associated color channel. Given the color channel 
histograms, the lower (bottom) and upper (top) color levels for each of the R, G, 
and B channels, denoted by Br, B a B B and Tr, T G , T B , respectively, are computed. 
It is noted these are not the actual highest and lowest values as outlying values 
may exist in the image being corrected. Rather, the lower and upper color levels 
for each channel represent the lowest and highest levels that are still consistent 
with the overall distribution of color levels in the image. More particularly, Br, B g , 
and B b , are each determined using the following equation: 

B-l B 

£"/<P#and £w, >pW, (20) 

/=o /=0 

where / is the color level in a range [0, L], N is the total number of pixels, and p is 
a small number chosen to eliminate the impact of noise at the top and bottom of 
the color level range. Similarly, T& T g , and T b are each determined using the 
following equation: 
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i-1 i-1 



and 



(21) 



i=t /=r+i 



In tested embodiments of the present system and process employing a pixel 
color level range of 0 to 255 integer levels, p was empirically chosen to be about 
0.5%. This represents an amount just slightly larger than one intensity level. It is 
noted that if more than 256 intensity levels are employed, the value of p should 
remain about 0.5%. However, less than 256 intensity levels are employed, then 
the value of/? should be adjusted to represent about 1 intensity level. For 
example, if 64 intensity levels are employed, the value of p should be about 
0.2%. 

Next, a histogram of the Y channel of the pixels in the image is computed, 
denoted by hy(r k ). The Y channel of a pixel, which is also called the luminance 
channel, is computed from the R, G, B color channels as follows: 



for an image having 256 levels in each color channel. For images having more 
or less levels, Y is computed proportionally. 

In addition to the Y channel histogram, the average color associated with 
each intensity Y level r k , denoted by c t =(R*, G*, B k ), is also computed. More 
particularly, computing c* for a particular intensity V level is accomplished as 
follows. First, the number of pixels in the image that exhibit an R channel having 
the intensity V level under consideration is divided by the total number of pixels 
in the image exhibiting that intensity Y level in any of the color channels to 
compute R*. Similarly, the number of pixels in the image that exhibit a G channel 
having the intensity V level under consideration is divided by the total number of 



76i? + 150G + 30S 
256 



(22) 
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20 



pixels in the image exhibiting that intensity /level in any of the color channels to 
compute G* . And finally, the number of pixels in the image that exhibit a B 
channel having the intensity V level under consideration is divided by the total 
number of pixels in the image exhibiting that intensity V level in any of the color 
channels to compute B* 

Given the Y channel histogram hy(r k ) and the average color associated 
with each intensity Y level c*, the average color of brightest pixels is computed as 
follows: 



L-l L-\ 

c ,op = (R,o P , G top , B lop ) = £ hy (r k )c k I £ h Y (r k ) 

*=/> k=p ' 



(23) 



where P is determined such that, 



L-l L-l 

15 2>/ <pN and H n i^P N 

i=p-\ i=p 



(24) 



Here again, in tested embodiments of the present system and process 
employing a pixel color level range of 0 to 255 integer levels, p was empirically 
chosen to be about 0.5%. 



An assumption is made in the design of the present system and process 
that the brightest pixels correspond to white/grey or glossy areas, and thus their 
color reveals the actual color of the light source. This is why the average color of 
the brightest pixels is computed. Further it is assumed that the overall image 
25 color is consistent with the color of the light source. Given this latter assumption, 
the R and B pixel color channels can be balanced with the G channel to improve 
the color of the image. The G channel is used as the reference for balancing 
because in a typical Bayer CCD chip, there are twice as many G pixels than R 
and B pixels, and the G value is considered to be more reliable. The first phase 
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of this balancing involves ensuring that the assumptions are valid, and to 
determine what, if any, balancing is required. To this end, it is first checked 
whether the R channel is balanced with the G channel. This involves the 
following considerations: 

1 • If \R, op - G lop \ < x (R lop + G lop )l2, or \B R - B G \ < x (B R + B G )/2, then the R and 
G channels are balanced and no change is required. In tested embodiments of 
the present system and process, x was set equal to 2.5%. This value was 
chosen empirically to provide a small tolerance range which even given noise 
there is a high confidence that the R and G channels are balanced. As such, 
other values for t should be chosen if the noise characteristics of the images 
dictate a larger value or allow for a smaller value. It is also noted that the value 
of t is dependent of the number of color levels, which in the case of the tested 
embodiments was 256. If more or less color levels are employed, the value of x 
should be adjusted proportionally. 

2. If R top > G top and B R > B G , then the dynamic range of the G channel is to be 
increased to match the R channel. 

3. If R top < G top and B R < B c , then the dynamic range of the R channel is to be 
increased to match the G channel. 

4. If none of the foregoing conditions apply, then the R and G channels are 
considered balanced because the image color is not consistent with the color of 
the brightest pixels. As the foregoing assumptions are not valid, no corrective 
action is taken. 

It is next checked whether the B channel is balanced with the G channel. 
This involves the following considerations: 
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1 . If \B lop - G top \ < x (B top + G lop )/2, or \B B - B G \ < x (B B + B G )/2, then the B and 
G channels are balanced and no change is required. 

2. If B top > G, op and B B > B c , then the dynamic range of the G channel is to be 
5 increased to match the B channel. 

3. If B top < G, op and B B < B G , then the dynamic range of the B channel is to be 
increased to match the G channel. 

10 4. If none of the foregoing conditions apply, then the B and G channels are 
considered balanced because the image color is not consistent with the color of 
the brightest pixels. As the foregoing assumptions are not valid, no corrective 
action is taken. 

15 It is noted that while in the foregoing procedure the R channel/G channel 

balance was assessed first followed by the B channel/G channel balance. This 
need not be the case. Rather the B channel/G channel balance could be 
assessed first followed by the R channel/G channel balance. In addition, the 
assumption that the G channel data is more reliable need not be adopted. In 

20 such a case, another of the color channels could be employed as the reference 
channel instead. 

If it is determined in the foregoing analysis that a color channel's dynamic 
range is to be increased to match another's range, this can generally be 

25 accomplished by computing the appropriate new top and bottom levels for the 

channel that is to be changed, and then linearly stretching the color range of that 
channel. More particularly, the following procedure is employed to stretch the 
color range of a channel. For the purposes of this description it will be explained 
how the B channel is increased to match the G channel. However, it is to be 

30 understood that the same procedure can be used stretch the dynamic color 
range of any channel to match another channel. 
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Given the previously determined values for by B G , T G and B B , T B : 



(l-ct)fl G and 



(25) 



5 



a(I-l) + (l-a)r 0 . 



(26) 



Eqs. (25) and (26) provide a way to control how aggressive the correction is by 
using a safeness parameter a to adjust the bottom and top levels of the color 

10 channel that is to be matched. No matching will be performed if a = 0, and the 
dynamic range of the color channel will be fully matched if a = 1 . This factor can 
be prescribed, or user specified. If prescribed, since this is an automatic tool, 
the value should be chosen so as to be conservative in the correction. This may 
not achieve the best correction, but should have the effect of not making any 

15 image look worse. For example, a prescribed value for a of about 0.75 would be 
a good conservative choice. However, a user specified value may be a better 
way go since there will be exceptions to the assumption that an image should fill 
the entire dynamic intensity range available because the photographer may want 
to achieve certain effects. Thus, a user specified safeness parameter provides a 

2 0 way to tailor the amount of correction to the circumstances of the image. 

Once the safeness parameter has been applied, a linear color range 
correction transformation is established and each old color level / of the B 

channel is mapped to a new color level f . Specifically, the transformation takes 
25 the form of: 



/ =al + b 



(27) 



where, 
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(t g -b g )/(t b -b b ) 

(b g T b -T g B b )/(T b - 
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